CDW, Superconductivity and Anomalous Metallic Behavior in 2D Transition Metal 

Dichalcogenides 
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We propose a theory for quasi-two-dimensional transition metal dichalcogenides that provides 
a unified microscopic picture of the charge density wave (CDW) and superconducting phases. We 
show, based on the electron-phonon coupling and Fermi surface topology, that a CDW order pa- 
rameter with six-fold symmetry and nodes (f-wave) gives a consistent description of the available 
experimental data. The elementary excitations in the CDW phase are Dirac electrons. The su- 
perconducting state has its origin on the attractive interaction mediated by phonons. The theory 
predicts strong deviations from Fermi liquid theory in the CDW phase. 



PACS numbers: 74.70.A, 71.45. Lr, 71.10.-w 

The quasi-two-dimensional (2D) transition metal 
dichalcogenides (TMD), 2H-TaSe 2 , 2H-TaS 2 , 2H-NbSe 2 
and 2H-NbS 2 have been extensively studied since their 
discovery jlj]. These are layered systems with a phase 
diagram where CDW order co-exists with superconduc- 
tivity at low temperatures (see Fig.[j]). Many physical 
properties such as resistivity, thermal expansion, sound 
velocity and so on, are very anisotropic and differ by or- 
ders of magnitude if measured parallel or perpendicular 
to the metallic planes ||||]. Neutron scattering studies 
have shown that these systems undergo a triple incom- 
mensurate CDW transition at a temperature T c d w with 
wave- vectors Qj = (1 — <5)bi/3 where are the three re- 
ciprocal lattice vectors for a triangular lattice with lattice 
spacing a (|bj| = 47r/(\/3a)) and 5 <C 1 is the incommen- 
surability (TaSe 2 also has a transition into a commensu- 
rate state, (5 = 0) j§. 

Although the phenomenological description of the 
CDW transition is well understood due to the work of 
McMillan M the current understanding of the micro- 
scopic origin of the CDW state is not complete. Fermi 
surface nesting, which is the main mechanism for CDW 
formation in ID systems (such as lT-TaSe 2 B) leads to a 
lattice deformation with momentum 2kp (where kp is the 
Fermi momentum), to the gapping of the Fermi surface 
and therefore to insulating behavior. Band structure cal- 
culations (BSC) do not show, however, strong nest- 
ing at the Fermi surface of TMD. Moreover, TMD be- 
come better metals (with a resistivity drop ||) inside of 
the CDW phase and the specific heat behaves like a power 
of the temperature, T, (instead of exponentially as in the 
case of an ID CDW transition). What is more striking 
is that TMD become superconductors at low tempera- 
tures. An alternative scenario was proposed by Rice and 
Scott (RS) based on earlier BSC M in which elec- 
tron scattering with momentum Qsp between different 
saddle-points (which produce a logarithmically divergent 
density of states) would drive the system into the CDW 
phase. The advantages of this picture are clear: it does 
not require nesting, leads to a gap at the saddle points 



and to the reduction of the resistivity. More refined BSC 
have shown that the saddle points are not as close to the 
Fermi surface as initially believed Q (see Figj^). 

With the advent of high temperature superconduc- 
tors (HTC) many experimental techniques have been 
perfected. Among them, angle resolved photoemission 
(ARPES) has nowadays momentum and energy resolu- 
tion that was unavailable when TMD were discovered. 
ARPES measures the electronic spectral function, the 
loci of the Fermi surface, and the momentum, k, and fre- 
quency, v, dependent electron self-energy, S(k, co). The 
imaginary part of S(k, ui) provides a direct measurement 
of the electron lifetime and for a Fermi liquid behaves 
like ^s[E(kp,uj — > 0)] cx 1/tq + boj 2 /Ep, where tq is the 
impurity scattering lifetime, Ep is the Fermi energy and 
b is a constant (we use units such that K = Hb = 1)- 
Recent measurements of S(k, oj) for HTC || have shown 
that S[X(fcp,aj)] tx |w| in clear disagreement with Fermi 
liquid theory but in agreement with the marginal Fermi 
liquid description |10(| . Stimulated by the recent devel- 
opments in ARPES, TMD have been extensively studied 
in the last two years. As a result of these studies a con- 
tradictory picture of the origin of the CDW state has 
emerged. Although there are indications for the opening 
of a CDW gap along the saddle point direction (but in no 
other place along the Fermi surface) pH , the measured 
|Qsp| is larger than 2k F and |Q,| (see Fig.|) p[. Fur- 
thermore, ARPES measurements of the electronic self- 
energy have shown that, analogously to the case of HTC, 
Q[S(fcp, lu — > 0)] oc l/r + 6|o;| (r|. These measurements 
clearly show that the current theories for CDW formation 
cannot describe the current experimental situation. 

We propose that, due to strong variations in the 
electron-phonon coupling and to imperfect nesting of the 
Fermi surface, the CDW gap is six-fold symmetric and 
has nodes (f-wave). This f-wave state has lobes (the 
largest value of the CDW gap, Max[|A c du,(k)|]) along 
the saddle point direction. This effect leads to the reduc- 
tion of the electron-phonon scattering in the same way 
as the RS mechanism. The low-lying excitations in the 
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CDW state are Dirac electrons associate with the nesting 
between different branches (around the T and K points) 
of the Fermi surface which are connected by Qi. Thus, 
the system remains a good metal in the CDW phase. 
Due to the loss of lattice inversion symmetry in the CDW 
phase, the Dirac electrons are coupled via a piezo-electric 
coupling to acoustic phonons. This coupling produces a 
damping of the Dirac electrons and leads to a S(k, oj) of 
the form (T <C T cdw ): 



S[£(fc F ,w^0)] = T^ 1 +aT + 1 \uj\ 



(1) 



for uj < Max[\ A c d w \] and 3[£(£;f,^ — * 0)] = 1/n, con- 
stant, for Maa;[|A cc j„,|] < u. We show that the data is 
in good agreement with these results. Finally, we argue 
that phonons drive the system to a superconducting state 
via a Kosterlitz-Thouless (KT) phase transition. We pro- 
pose that critical fluctuations of the superconducting or- 
der parameter, A s , lead to a pseudo-gap behavior above 
the transition temperature, T c , to a drop of the magnetic 
susceptibility, and to diamagnetic response above T c (as 
seen experimentally fTi]), 

Unlike ordinary metals, the electron-phonon coupling, 
g\ (k, q) , (where q = k' — k is the momentum transfer 
for an electron being scattered by a phonon from k to 
k' and A is the phonon polarization) in transition metals 
is highly dependent on the electron momentum. This 
dependence is responsible for many anomalies observed 
in the phonon spectrum in these systems [fj"5f . Numerical 
studies have shown that g\(k., q) can vanish at certain 
points in the Brillouin zone fl6|| . A simple tight-binding 
calculation gives, 



5A(k,q) 



V it. 

x cos[(k + q/2) • n 



^C|„|e A:q • nsin(q- n/2) 



(2) 



where N is the number of atoms, M is the atom mass, 
^A,q is the phonon frequency, e\ q the polarization vec- 
tor (e^ _ q = e>, q ), n the nearest neighbor vector and C| n | 
depends only on |n|. It is easy to see that g\(k, Qi) van- 
ishes along the nodal lines in Fig.||. BSC (and ARPES 
measurements) show that the Fermi surface is hole-like, 
centered around the T and K points, and intercepts 
the nodal lines at the Dirac points (see FigJ||). At 
these points the CDW order parameter vanishes since 
Acdiu(k) <x g\(k,Qi). Observe that Qi connects nodal 
points located at different branches of the Fermi surface. 
Moreover, as shown by numerical calculations || these 
points are well nested and a CDW state can emerge. In a 
triple CDW phase, however, the situation is more compli- 
cated because Q2 connects the nodal points in the same 
branch of the Fermi surface. This would imply a finite 
coupling at these points and the opening of a gap there 
as well. We observe, however, that in this direction the 
Fermi surface is not nested [pi in contrast with the case 



discussed above (these details are not shown in FigJ|). 
Thus, a gap cannot open at the Dirac points. The con- 
clusion is that the CDW gap will change sign in six points 
along the Fermi surface leading to the situation depicted 
on Fig.g. This is the f-wave CDW state. 

Let Ck,cr (c t ) be electron annihilation (creation) op- 
erator for an electron with momentum k and spin pro- 
jection a (a =t, I). We can define the spinor operators 

*L(k) = (4,^4 +Qi , CT ) = (V>W(k),< ilCT (k)) (3) 

where +, — indicates if the fermion is particle or hole 
(anti-particle) like. The Hamiltonian for the excitations 
close to the nodal points is 

H d=Y. *U k ) Wk^i * + v o k \\,i aX ] **.*( k ) ( 4 ) 



where vf is the Fermi velocity, vq = |9fcA c d UJ (k)|, k± t i 
(fen i) is the momentum perpendicular (parallel) to the 
Fermi surface . All these quantities are calculated at the 
Dirac points. Here a a with a — x,y,z are Pauli ma- 
trices which act in the spinor subspace. Hamiltonian 
([I|) describes a system of Dirac fcrmions with energy 

e^ k = ±ej, k where e;, k = ^ v F k l,i + u o fc {[. 4 - In wnat 
follows we drop this index i since the Dirac electrons are 
decoupled. 

The triple CDW transition breaks the inversion sym- 
metry of the lattice leading to a three sublattice structure 
jl| . In this case a piezo-electric coupling can develop be- 
tween the Dirac fermions and acoustic phonons (ItJ . The 
acoustic phonons are described by 



H A = ^-jdr [n 2 (r) + plcl{VHv) f] 



(5) 



where c s is the sound velocity, pl is the lattice mass den- 
sity and II(r) and $(r) are the momentum and phonon 
fields, respectively. The electron-phonon coupling is 



Ho =K 



dr*(r)53*t( r )¥ a ( r ) 



(6) 



where k is the piezo-electric coupling constant. The 
electron Green's function for the problem described by 
Hamiltonian (§, (g) and (§) can be calculated self- 
consistently flSfl . The main result of the interaction is the 
damping of these excitations. As a result 3[E(fci?,w)] is 
given in (@) with t^ 1 = 27 meV, a = 2.14 and 7 = 0.212. 
These results are valid for li <C Max [|A ce fa,|] ~ 60 meV 
(obtained from the onset of the optical absorption 
and ARPES For uj > Max[\A cdw \] the Dirac elec- 

tron description breaks down and Sy[S(A;ir, u)] ~ t^ 1 = 
73 meV. In Figj3| we plot our results against the exper- 
imental data. The agreement is remarkable given the 
simplicity of the model. Furthermore, the predictions of 
the marginal Fermi liquid phenomenology originally pro- 
posed for the HTC [JlOt] can now be extended to the case 
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of MTD. In fact, recent dynamical transport measure- 
ments have pointed out the striking similarities between 
the MTD and HTC @. 

Besides leading to damping, the phonons generate a 
retarded interaction. This retarded interaction, like in 
an ordinary superconductor, leads to pairing in the sin- 
glet channel. After tracing out the phonons the pairing 
Hamiltonian becomes: 

H P = -g ]T < 6 <„< T (-k')VcT (k)Vda (-k) (7) 

k,k' 

where g ~ k 2 /ujjj is the coupling constant [ojjj « c s A, 
with Aw 1/a, is a Debye frequency). At the mean- field 
level the pairing Hamiltonian is written as 



Hp = 



k,a,6 



< 6 A,^ )T (k)^ 4 (-k) + /i.c. 



where 



A, 



E 



(8) 



(9) 



is the superconducting order parameter. The problem 

described by (0) and (|J) reduces to the diagonalization 

of a 4 x 4 matrix via a Bogoliubov transformation. The 
?± - j- 1,2 , i A p. in tne 



eigen-energies are: 



El 



superconducting phase the Fermi surface is fully gapped 
by a CDW gap along the T-K direction and a supercon- 
ducting gap along the T-M direction. |A S | given by: 



\A s (T,g)\ = 2Tcosh~ 1 



cosh I — — 1 e T s 



(10) 



V T 9c 

where g c = Att 3 ^ 2 ^/vfVo/A. At T — we have: 

\A s (0,g)\=iTrv F v (l/g c -l/g) (11) 

which shows that superconductivity is only possible for 
g > g c at T — 0. Thus, g — g c is a quantum critical point 
(QCP) and g c is the critical coupling constant. This re- 
sult implies that there is a critical lattice spacing, a c , be- 
low which superconductivity is not possible. The critical 
mean-field temperature is T*(g) = |A S (0, g)\/{2 ln(2)). 
Notice that the dependence of the order parameter with 
the coupling constant is very different from the BCS ex- 
pression (which does not require a critical coupling con- 
stant). This transition, as was noted earlier [|Fj), can be 
associated with the spontaneous chiral symmetry break- 
ing in the Yukawa-Higgs model. In our case, however, 
the order parameter is complex implying that there is an 
extra U(l) symmetry associated with the phase of the 
order parameter. Because of the strong phase fluctua- 
tions in 2D the phase transition can only be of the KT 
type (2D-XY) due to the unbinding of vortex-anti- vortex 
pairs at a temperature 



TKr(g) = na s (g)/{2m* 



(12) 



where m* is the pair effective mass and a s (g) cx 
|A S (0,<?)| 2 is the planar superfluid density. This tran- 
sition, however, does not produce true long-range or- 
der. The weak coupling between layers changes the uni- 
versality class of the transition to 3D-XY with T c w 
T KT + b/\n 2 {T KT /(cU ± )) where b is a number of or- 
der unit, U± is the coupling energy per unit of length 
between layers, and c is the inter-layer distance (we as- 
sume, c <C Tkt/U±). Thus, T c grows with decreasing c. 
As we can see the relevant parameter which controls the 
superconducting transition is a/c. In Fig.l we plot the 
transition temperatures as a function of this parameter. 
We can clearly see the anti-correlation between T c d w and 
T c . Notice that TaSe 2 with T c w 0.1 K is close to the 
QCP. The growth of T c under pressure J2l[] is also con- 
sistent with our picture. Because T* > T c we expect the 
a pseudo-gap region for T c < T < T* where the order 
parameter is developed but there is no phase coherence. 

As we have shown there are many similarities between 
the problem discussed here and HTC. Perhaps one of the 
most interesting is the anti-correlation between the T c d w 
and T c . This is the same type of anti-correlation observed 
between the pseudo-gap energy scale and T c in the case 
of HTC p2| . The main difference, of course, is that at 
T c dw there is a true second order phase transition while 
the pseudo-gap temperature seems to be a crossover en- 
ergy scale. In fact, the metallic state described here is 
similar in many respects to the nodal liquid description 
of HTC (2^]. However, the MTD Dirac fermions are not 
related superconductivity but to a CDW state. In TMD 
the anomalous metallic behavior is due to the coupling 
to acoustic phonons (which are critical modes associated 
with the breaking of translation symmetry) while in the 
case of HTC this is not so clear. As shown in ref. [ p"8[ , 
phase fluctuations alone are not capable of producing 
marginal Fermi liquid behavior. In order to obtain this 
kind of behavior the system needs be close to a different 
type of QCP. There are many attempts to generate hid- 
den QCP for HTC. Among them, the idea of stripes (a 
state closely related to CDW) have been discussed in the 
literature [|4|. It is worth noticing that TMD also show 
stripe phases associated with the breaking of the hexag- 
onal symmetry in the triple CDW phase pa]. Another 
possibility is related with the breaking of the time rever- 
sal symmetry and generation of gapped superconducting 
phases . More recently it has been proposed that due 
to strong electron-electron correlations an exotic CDW 
state might be the hidden QCP of HTC (2(| and that 
the pseudo-gap temperature marks a real second order 
transition that is smeared by disorder. If this is indeed 
correct then the anomalies in these 2D systems can have 
the same origin. 

In summary, we have proposed a microscopic theory 
for the CDW and superconducting phases of TMD that 
involve the formation of a gapless CDW state. The ele- 
mentary excitations are Dirac fermions that pair up due 
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to the coupling to phonons and generate a superconduct- 
ing state at low temperatures. We show that in the CDW 
state the coupling to phonons leads to marginal Fermi 
liquid behavior and to many anomalies in the physical 
properties. We show that our theory can explain the 
available experimental data including ARPES. We pre- 
dict that high-resolution ARPES should be able to mea- 
sure the modulation of the CDW gap around the Fermi 
surface and the opening of a superconducting gap at the 
Dirac points. The similarities between the phenomena 
described here and the one observed in other layered sys- 
tems (such as HTC) is striking and the study of such ex- 
otic CDW might illuminate the way for the understand- 
ing of the complex physics in transition metal systems. 
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FIG. 1. Phase diagram: Stars, T c ^ w ; Filled squares, T c . 
From left to right, TaSe 2 , TaS 2 , NbSe 2 and NbS 2 |J|. a is 
the in-plane lattice spacing and c is the inter-layer spacing. 




FIG. 2. Skematic plot of the Fermi surface accordingly to 
||. Dashed lines: nodal lines associated with Qi; Filled cir- 
cles: Dirac points; Empty circles: Saddle points; Thick line: 
proposed CDW gap. 
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FIG. 3. 9[E(fc F ,w)L Experiments: squares (T = 34K), 
pentagons (T = 76K) Q; Theory: dotted (T = 0), continous 
(T = 34K), dashed (T = 76K). 
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